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ABSTRACT 

We detect the likely optical counterpart to previously reported X-ray QPOs in 
spectrophotometry of the intermediate polar GK Per during the 1996 dwarf nova out- 
bm'st. The characteristic timescales range between 4000-6000 s. Although the QPOs 
are an order of magnitude longer than those detected in the other dwarf novae we 
show that a new QPO model is not required to explain the long timescale observed. 
Wc demonstrate that the observations are consistent with oscillations being the re- 
sult of normal-timescale QPOs beating with the spin period of the white dwarf. We 
determine the spectral class of the companion to be consistent with its quiescent clas- 
sification and find no significant evidence for irradiation over its inner face. We detect 
the white dwarf spin period in line fluxes, V/R ratios and Doppler-broadened emission 
profiles. 
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1 INTRODUCTION 

' GK Per (Nova Per 1901; CampbeU 1903), belongs to a sub- 
group of cataclysmic variables (CVs) called Intermediate 
Polars (IPs). In these systems an asynchronously-rotating, 
, magnetic white dwarf accretes material from a less-massive, 
' late-type companion filling its Roche lobe. Gas leaving the 
companion star attempts to form an accretion disc around 
the primary star but its magnetic field either prevents the 
formation of the disc or truncates it near the white dwarf. 

GK Per was identified with the X-ray source A0327-I-43 
by King, Ricketts & Warwick (1979) and confirmed as an IP 
by the detection of a 351 s X-ray spin pulse by Watson, King 
& Osborne (1985; hereafter WKO) and Norton, Watson & 
King (1988). The same period was subsequently found in 
optical photometry by Patterson (1991). GK Per has the 
longest orbital period from the sample of known CVs, Pqj.]-, 
= 2d, (Crampton, Cowley & Fisher 1986; hereafter CCF). 
The wide binary separation combined with a relatively weak 
magnetic field (~ 1 MG) means that a truncated accretion 
disc must be present if current theories of disc formation 
are correct (Hameury, King & Lasota 1986). The presence 
of a disc has yet to be confirmed by direct observation, al- 
though the system does undergo dwarf nova outbursts every 
2-3 years where its optical brightness increases from 13th to 
10th magnitude (Sabbadin & Bianchini 1983). The most- 



likely mechanism for dwarf nova outbursts is a thermal in- 
stability within an accretion disc (Osaki 1974) . GK Per out- 
bursts have been modelled as such by Cannizzo & Kenyon 
(1986) and Kim, Wheeler & Mineshige (1992). 

This paper is a continuation of paper I (Morales- Rueda, 
Still & Roche 1996), in which we presented spectrophoto- 
metric observations of GK Per taken on the rise to its 1996 
outburst (Mattel et al. 1996). We reported the detection 
of quasi-periodic oscillations (QPOs) within the Doppler- 
broadened emission lines of H I and He ll. This provides 
an opportunity to map the velocity structure of the oscilla- 
tions. QPOs are defined as low-coherence brightness oscil- 
lations thought to be associated with material within the 
inner accretion flows of CVs. Theoretical models developed 
to explain QPOs consider the presence of dense blobs of ma- 
terial orbiting in the inner regions of the accretion disc (Bath 
1973), or non-radial pulsations over the surface of the white 
dwarf (Papaloizou & Pringle 1978), or radially-oscillating 
acoustic waves in the inner disc (Okuda et al. 1992; Godon 
1995). In these models, the QPO timescales match observa- 
tions of dwarf novae and are of the order of a few hundred 
seconds. However the QPO periods detected in GK Per are 
an order of magnitude longer than this. Previous to this 
paper they have only been detected in X-ray data taken 
during outbursts; WKO discovered them in 1.5-8.5 keV EX- 
OSAT data at the peak of the 1983 outburst, while Ishida 
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Figure 1. The visual light curve of GK Per during the 
; 996 outburst obtained from the Var iable Star Network 
(http://www.kusastro.kyoto-u.ac.jp/vsnet/). The arrows indi- 



cate the times of our spectrophotometric observations. 



et al. (1996) report a second detection at 0.7-10 keV with 
ASCA during the rise to the 1996 outburst discussed in this 
paper. 

To explain the long timescales WKO suggested the 
QPO mechanism is caused by beating between the 351 s 
white dwarf spin period and inhomogeneous gas orbiting at 
the inner edge of the accretion disc. HeUier & Livio (1994; 
hereafter HL) noted that the X-ray hardness ratio varies 
over the QPO cycle as expected from photoelectric absorp- 
tion by cool gas and that a period of a few thousand sec- 
onds is consistent with the orbital frequency of gas if it is 
deposited onto the disc by a gas stream which has partially 
avoided impacting the outer disc rim and follows a ballistic 
trajectory. They propose that the QPO mechanism is X- 
ray absorption by vertically-extended blobs of gas orbiting 
at this preferred inner impact radius. In paper I we deter- 
mined that the characteristic velocity structure of the optical 
counterpart to the QPOs observed by Ishida et al. (1996) is 
consistent with blobs in the inner disc. In the current paper 
we present further analysis which indicates that the optical 
QPO is also driven by absorption, but favours strongly a 
beat model over the disc-overflow interpretation. 



2 OBSERVATIONS 

Between 1996 February 26 and 28, 6-8 days before the 
ASCA pointings of Ishida et al. (1996), we obtained spec- 
trophotometry of GK Per using the Intermediate Dispersion 
Spectrograph mounted on the 2.5 m Isaac Newton Telescope 
(INT) on La Palma. Table 1 gives a journal of observations. 
In Fig. 1 we show a visual light curve obtained by the Vari- 
able Star Network during the 1996 outburst, with arrows 
indicating the days on which we made observations. The 
quick readout mode was used in conjunction with a Tek- 
tronix CCD windowed to 1024 x 150 pixels to reduce dead 



Table 1. Journal of observations. E is the cycle number plus 
binary phase with respect to the ephemeris given by Crampton, 
Cowley & Fisher (1986). Phases have been adjusted by tt/2 so 
that phase corresponds to superior conjunction of the white 
dwarf. 



Date Start End 

(1996 Feb) (UT) 



Start End No. of 

{E- 2 000) spectra 



26 
27 
28 



20.06 
20.04 
20.06 



0.12 
0.07 
0.04 



617.119 
617.620 
618.121 



617.204 
617.704 
618.204 



109 
116 
117 



time and obtain good sampling of the spin cycle. The expo- 
sure times and resolution of the data were already described 
in paper I. 

After debiasing and flat-fielding the frames by tungsten 
lamp exposures, spectral extraction proceeded according to 
the optimal algorithm of Home (1986). The data were wave- 
length calibrated using a CuAr arc lamp and corrected for 
instrumental response and extinction using the flux stan- 
dard HZ 15 (Stone 1977). The spectrograph slit orientation 
of PA 249.1° allowed a 15th magnitude nearby star approx- 
imately 0.5 arcsec ENE of GK Per to be employed as cali- 
bration for hght losses on the slit. 

We also have available to us spectroscopy of various K- 
type stars from 1995 October 11 to 13 obtained from the 
INT and from 1995 October 30 to November 2 with the 
2.1 m telescope in the McDonald Observatory in Texas. The 
INT instrumental setup was identical to the one used for 
the 1996 observations described above. For the McDonald 
data, the low-to-moderate resolution spectrometer ES2 was 
employed in conjunction with the Til CCD and a grating 
ruled at 1200 lines mm~^ covering the wavelength region 
A4196 A-A4894 A giving a resolution of 200 km s"^ at H/3. 

The spectra were flat-fielded, optimally extracted and 
wavelength calibrated also in the standard manner. Flux 
calibrations were applied using observations of the standards 
HD19445 (Oke & Gunn 1983) and Feige 110 (Stone 1977) for 
the October and November data respectively. Table 2 gives 
a list of the K-type templates observed over both runs. 



3 RESULTS 

3.1 Average spectra 

Fig. 2 presents the average of all the data collected on 1996 
Feb 28. It is characterised by a flat continuum, broad Balmer 
and He I lines in emission, high excitation lines of Hell, Nlll 
and Cm and numerous faint, narrow absorption features of 
Fe I, Cal, Tin and Sril that had been identified as signa- 
tures of the K-type secondary star by Kraft (1964), Gal- 
lagher & Oinas (1974), CCF and Reinsch (1994). 

We employed K star spectral templates to determine 
which luminosity class best matched the secondary star in 
this system during outburst and search for signatures of in- 
creased X-ray irradiation. Using CCF's fit to the orbital ra- 
dial velocity of the secondary star we shifted out the orbital 
motion of the absorption lines with a quadratic rebinning 
algorithm. We binned in velocity the spectra of GK Per and 
the K-type templates to ensure that they all had identical 
wavelength ranges and dispersions. We employed the op- 
timal subtraction algorithm of Marsh, Robinson & Wood 
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Figure 2. The top spectrum is an average of all GK Per spectra obtained during the third night of observations. The bottom spectrum 
corresponds to the KlIV template HD197964 multiplied by 3x10"*. The middle spectrum is the residual resulting from the subtraction 
of the template from the averaged GK Per spectrum and probably resembles the spectrum of the accretion flow. 



(1994) to determine the K star spectral type - we multiply 
the template by a monochromatic constant which represents 
the contribution to the spectrum from non-stellar sources of 
light and subtract the resulting spectrum from the GK Per 
data. The residual was smoothed using a high-pass band fil- 
ter (FWHM of gaussian = 13 A), and a test performed 
between the original and smoothed residual. This is an it- 
erative procedure to determine the optimum value of the 
monochromatic constant which continues until is min- 
imised. Table 2 lists the templates, their spectral classes, 
and the reduced obtained after applying optimal sub- 
traction. The best fit template is the Kliv star HD197964 
which provided a reduced of 2.5. The secondary star con- 
tributes 13 per cent of the total hght in this spectral region 
on the third night of observations. This compares to 33 per 
cent found by CCF and Gallagher & Oinas (1974) during 
quiescence indicating that the accretion flow has increased 
in brightness. 

The best-fit luminosity classes are consistent with the 
quiescent classifications of K2ivp by Kraft (1964), K2ivp- 
K2v by Gallagher & Oinas (1974), KOlll by CCF, and K3v 
by Reinsch (1994). We find the spectral type to be con- 
stant across our two phase samples - one during which a 
large area of the white dwarf-facing surface is visible and the 
other when it is mostly limb-occulted. Consequently there 



Table 2. A list of template K star used to determine the best-fit 
spectral type for the secondary star. 



Name 


Spectral 




Name 


Spectral 






type 






type 




October 1995 




November 1995 




HR190 


Kim 


4.3 


13 Lac 


KOni 


4.0 


HR8688 


Kliii 


6.6 


1 Peg 


Klni 


5.7 


HR8415 


K2in 


6.5 


69 Aql 


K2ni 


3.9 


HR8632 


K3in 


6.2 


39 Cyg 


K3ni 


4.0 


HR8974 


Kliv 


5.9 


33 Vul 


K3.5III 


3.8 


HR8881 


Klv 


6.2 


3r)Cep 


KOiv 


5.9 


HR222 


K2v 


3.0 


HD197964 


Kliv 


2.5 


HR8832 


K3v 


3.2 









is no observational evidence for an increase in irradiating 
flux from the accretion regions over the inner face of the 
companion star, although we are limited by a small range of 
spectral templates and poor orbital sampling. 

In order to measure integrated emission line fluxes 
from each of the three nights, we fitted a third order 
polynomial through wavelength bands relatively free of 
line features (AA4147-4212A, AA4278-4306A, AA4560-4608A, 
AA4770-4838A and subtracted the fit from the data. Fluxes 
were measured by summing under each line profile and these 
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Table 3. Emission line fluxes in units of 10~^^ erg cm~^ s~^ from the spectra of GK Per on 1996 Feb 26-28. The minimum and maximum 
fluxes throughout the night as well as the average nightly flux are provided. The error on the average flux is the standard deviation. 
Negative fluxes correspond to line absorption. 

Line Flux Average Flux Average Flux Average 

range flux range flux range flux 



26/2/96 27/2/96 28/2/96 



H7 


1.24 - 


2.97 


2.05 


± 


0.34 


1.57 - 


3.31 


2.27 


± 


0.30 


2.09 


- 4.03 


3.06 


± 


0.37 


H/3 


1.43 - 


3.63 


2.32 


± 


0.40 


1.73 - 


- 3.27 


2.45 


± 


0.35 


2.17 


- 4.57 


3.41 


± 


0.46 


He I A4387.9A 


-0.70 - 


- 0.31 


-0.04 


± 


0.15 


-0.18 - 


- 0.29 


0.03 


± 


0.08 


-0.14 


- 0.35 


0.11 


± 


0.09 


He I A4437.6A 


-0.83 


- 036 


-0.29 


± 


0.18 


-0.60 - 


- -0.15 


-0.34 


± 


0.09 


-0.68 


- -0.23 


-0.44 


± 


0.09 


He I A4471.7A 


-0.47 - 


- 0.72 


0.13 


± 


0.18 


-0.36 - 


- 0.41 


0.12 


± 


0.14 


-0.28 


- 0.68 


0.25 


± 


0.18 


He I A4713.2A 


-0.70 - 


- 0.52 


0.02 


± 


0.14 


-0.18 - 


- 0.27 


0.05 


± 


0.08 


-0.22 


- 0.31 


0.02 


± 


0.10 


He I A4921.9A 


-0.06 - 


- 0.93 


0.38 


± 


0.15 


0.18 - 


- 0.58 


0.39 


± 


0.09 


0.26 


- 1.03 


0.68 


± 


0.13 


Hen A4541.7A 


-0.46 - 


- 0.25 


-0.06 


± 


0.14 


-0.30 - 


- 0.30 


0.02 


± 


0.10 


-0.21 


- 0.26 


0.02 


± 


0.11 


Hen A4685.8A 


4.01 - 


7.36 


5.23 


± 


0.77 


3.24 - 


- 7.72 


5.44 


± 


0.78 


6.04 - 


- 12.19 


8.27 


± 


0.94 


Bowcn blend 


0.04 - 


1.05 


0.50 


± 


0.17 


0.39 - 


- 1.08 


0.65 


± 


0..13 


0.71 


- 1.89 


1.06 


± 


0.17 



Table 4. Power indices a obtained from fitting with a power 
law function to the integrated flux of each emission line and the 
continuum over the three nights of observations. 



Line 



Power index Line 



Power index 



H7 

HeiA4388A 
He I A4472A 
He II A4686 A 
Continuum 



2.3 ± 0.1 
1.1 ± 0.4 
2.1 ± 0.7 
3.5 ± 0.1 
2.41 ± 0.01 



H/3 

He I A4438A 
He I A4922A 
Bowen blend 



2.7 ± 0.1 
2.5 ± 0.8 
3.5 ± 0.5 
2.3 ± 0.2 



are provided in Table 3. The intensity of the continuum, the 
lines and the relative intensity of He 11 A4686 A with respect 
to the Balmer lines, increases from the first night to the last 
as the system approaches the outburst maximum. We fit 
the emission lines during the three nights with a power law 
function of time F ^ t" , and provide the index a for each 
line in Table 4. 

Power-law fits of the form = on each consecutive 
night provide a = -1.61 ± 0.03, -1.43 ± 0.03 and -1.39 
± 0.11. Continuum slope changes slightly within statistical 
uncertainties during the observing run, the spectra becom- 
ing bluer with time consistent with a rise in temperature 
through the accretion flow. These indices are inconsistent 
with an accretion disc emitting as a discrete set of black- 
bodies (Pringle 1981). 

A comparison of the Feb 28 averaged spectrum with the 
spectra presented by Reinsch (1994) reveals that the Balmer 
line fluxes are ~1.7 times larger than during quiescence 
and the He 11 A4686 A feature and the C iii/N iii A A4640-50 A 
Bowen blend are ~5.3 times brighter. He I A4471.7A and 
He I A492I.9A are 1.4 and 2.3 times brighter during this 
outburst stage, respectively. In quiescence the Balmer lines 
are the brightest emission lines, whereas the strongest line 
in the current data is He 11 A4686 A. Szkody, Mattel & Mateo 
(1985) and CCF present spectra of GK Per taken during the 
1983 outburst maximum and 20 days after outburst respec- 
tively in which this behaviour is also clear. 



3.2 Radial velocities 

In paper I we provided an analysis of the emission line ve- 
locities. To complete the radial velocity analysis we now 



consider the absorption lines. In Sec. 3.1 we determined 
that our best secondary star template has a spectral type 
of Kliv. By masking out the emission lines in individual 
GK Per data and subtracting flts to the continua from all 
spectra, we were able to cross-correlate the absorption spec- 
trum of GK Per with our template (Tonry & Davis 1979). 
We corrected the resulting radial velocities by the systemic 
velocity of the template star (-6.5kms~^; Evans 1979) and 
fitted them with a circular function: 



V = 'y + K sin 27r [0 — 1 



(1) 



Orbital phases were adopted relative to the corrected CCF 
ephemeris, where <j)o corresponds to superior conjunction of 
the white dwarf. 7 represents the systemic velocity of the 
binary, K is the radial velocity semiamplitude of the com- 
panion star and (j> is the orbital phase. 

We combined the radial velocities measured by previ- 
ous authors (Kraft 1964; CCF; Reinsch 1994) with our own 
values and plot them together in Fig. 3. We assume that 
the errors on all individual measurements previous to this 
study are equal to the mean error of 20 kms~^. The solid 
curve is the fit to all data, providing 7 = 30 ± lkms~^, 
K = 119±2kms"^and <^o = 0.998 ±0.003. The dot-dashed 
curve is a fit to all the data excluding the current set where 
7 = 22±2kms-\ K = 128±2 kms'^and 0o = 0.009±0.003, 
providing reasonable agreement although the fits are not 
consistent within the given errors. Martin (1988) showed 
that an elliptical fit can account approximately for irradi- 
ation processes over the inner face of the secondary star. 
Elliptical fits to the quiescent data have already been pro- 
duced by CCF and Reinsch (1994). We do not have suitable 
phase sampling to produce a significant elliptical fit with 
the current data. Therefore although we find no evidence 
for secondary star irradiation in the absorption line radial 
velocities during outburst, our phase coverage prevents us 
from ruling it out. 



3.3 Emission line profiles 

The continuum-subtracted data are presented as time-series 
of selected line profiles over each night of observation in 
Fig. 4. At times theses lines display double-peaked profiles. 
This is often considered a signature of accretion disc emis- 
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Figure 3. Radial velocities obtained from the current data by cross-correlating absorption features from the companion star against a 
Kliv template, combined with similar measurements during quiescence from Kraft (1964), CCF, and Reinsch (1994). We provide two 
circular fits to the data. 



sion (Smak 1981), however the velocity structure across the 
accretion flow is so complex in IPs that this cannot be con- 
sidered a conclusive detection of the accretion disc. The pro- 
files are asymmetric where the peak apparently shifts from 
the blue to the red and back again in the Balmer lines over 
the observations. This behaviour is reminiscent of emission 
from a localised region in the system such as the bright spot 
where the accretion stream strikes the outer rim of the disc, 
or an irradiated region on the secondary star, but we note 
that the orbital phasing of the observed H/? feature is in- 
consistent with both interpretations. Moreover, the orbital 
phases at which we see these variations are not those at 
which the hot spot and the irradiated face of the secondary 
are best observed, i.e. phases 0.8 and 0.5 respectively. The 
profile variations of the Hel and He ll A4686 A lines are dif- 
ferent to those of H/3 either because they originate from 
different locations or are more sensitive to intervening ab- 
sorption regions. 

The most interesting variation occurs in the blue wings 
of all these profiles. First we note that the profiles are asym- 



metric about their rest velocities, regardless of orbital phase, 
where each line has a red bias. This shift is much larger than 
the systemic velocity of the binary, measured from secondary 
star photospheric lines. Secondly we note that this asymme- 
try is periodic, at least in H/3 and HellA4686A, and this 
period corresponds to the kilo-second QPOs we reported in 
paper I. The QPOs manifest in blue-shifted material and ap- 
pear to be the result of absorption either of the line source 
or the underlying continuum. We have presented trails of 
H/3 and HellA4686A against QPO phase after subtracting 
the nightly average from each spectrum in paper I. 

After the orbital period, the third likely signal present 
in these trails is the 351 s spin period of the white dwarf. We 
attempted to remove the orbital variations in the line pro- 
file by shifting out the motion of the white dwarf according 
to the ephemeris and radial velocity fit of CCF. The QPO 
contribution was accounted for approximately by combining 
the resulting spectra into 40 bins phased over the QPO cycle 
and subtracting the spectrum in the bin nearest in time from 
each individual spectrum. Fig. 5 shows the resulting trails 
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Figure 4. Trailed spectra of H/3, HeiiA4686A and tiie average of HeiA4388A, HeiA4472A, HeiA4713A, HeiA4922A). Each row- 
corresponds to a different night, the top row being the first night of observations. Images have separated linear greyscales where black 
corresponds to emission. 



of H/3, He II A4686 A and the sum of the Hel lines binned 
into 30 bins, over the spin period using Ishida's et al. (1992) 
ephemeris. The modulated signal is faint during the Feb 26 
and 27, but clearly present in the H/3 and He ll A4686 A pro- 
files on Feb 28 extending out as far as ~1000kms~^ in 
the H/3 profile. We see one modulation per cycle with signal 
moving from the red peak to the blue. This is reminiscent of 
the spin signal found in the trails of the IP RX J0558+5353 
(Still, Duck & Marsh 1998), although in that case most of 
the power occurred on the 1st harmonic of the spin period, 
indicating accretion onto two poles of the primary star. In 
the current trails of GK Per we see no evidence for power 
on the 1st harmonic. Harlaftis & Horne (1998) postulate 
that the origin of this spin pulsed emission is the region 
where the disc material is threaded onto the magnetic field. 
Similar trails but for spectra obtained during quiescence 
are presented by Reinsch (1994), where the spin signal is 
also clear on the fundamental frequency in the Balmer and 
HellA4686 A lines but not very strong in He I. 



3.4 V/R ratios 

By measuring the fluxes under the continuum-subtracted 
blue wings (from -1100 to kms~^) of the H/3 and 
He II A4686 A lines and dividing these by the flux under the 
red wings (from to 1100 kms~^) we produce a time-series 
of V/R ratios. These are plotted in Fig. 6 for the three nights 
of observation. Kilo-second oscillations are observed over all 
three nights. A power search over the ratios was performed 
using the Lomb-Scargle algorithm (Scargle 1982) and the 
QPO periods found are listed in Table 5. The errors quoted 
are only an estimate of the minimum error and depend on 
the frequency sampling. A significance test (a variant of the 
randomisation Monte Carlo technique, Linnell Nemec & Ne- 
mec 1985) was run by iteratively searching for periods after 
small shifts of the data had been performed. After 1000 per- 
mutations we found that the periods present in the V/R 
ratios are within the quoted errors with 95 per cent confi- 
dence. Note that the QPOs tend to shorter periods over the 
three nights. 

In order to determine the nature of higher-frequency 
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Figure 5. Trailed spectra of H/3, Hell A4686 A and the average of Hel (as in Fig. 3) binned in the spin period. Note that the spin cycle 
is repeated twice. As in Fig. 3, each row represents a different night, the top row being the first night. Each image is on a separate linear 
intensity scale. 



Table 5. QPO periods measured from the V/R ratios of H/3 and 
Hen A4686 A. 



1996 Feb 



QPO 



1996 Feb 



QPO 



H/3 26th 6362 ± 3 
H/3 27th 5025 ± 2 
H/3 28th 3964 ± 1 



He II A4686 A 26th 6568 ± 3 
He II A4686 A 27th 5241 ± 2 
He II A4686 A 28th 4333 ± 1 



variations we created a coarse version of each V/R curve 
over 40 time bins on each night, each bin approximating 
to one white dwarf spin cycle. The QPO signal was filtered 
out by subtracting the bin nearest in time from each V/R 
measurement. We searched for power in the modified V/R 
ratios using the Lomb-Scargle algorithm and obtained the 
power spectra plotted in Fig 7. 

Four of the spectra clearly show power at the white 
dwarf spin period (246 cycles day~^), but we cannot de- 
termine whether there is any power on the 1st harmonic 
which occurs beyond the Nyquist limit. We folded the mod- 
ified V/R ratio data using the spin ephemeris from Ishida 



et al. (1992) using 40 bins and plot them in Fig. 8. The V/R 
ratios show sinusoidal behaviour but not as clearly as in Gar- 
lick et al. (1994) and Reinsch's (1994) quiescence data. The 
maximum in these curves has previously been observed in 
quiescence at spin phase rather than phase 0.25. 



3.5 Emission line and continuum fluxes 

We conducted a power search across the same regions of 
continuum listed in Sec. 3.1 and the integrated flux over 
each emission line. Power is present at kilo-second periods 
whose values decrease on consecutive nights; see Table 6 and 
in Fig. 9 we present the power spectra for the integrated 
fluxes of H/9, He ll A4686 A and the continuum. No signifl- 
cant power was found in the spin period during any of the 
nights. Garlick et al. (1994) and Reinsch (1994) flnd clear 
modulations in the intensity of the Balmer lines with the 
spin period during quiescence but no signiflcant contribu- 
tions at the spin period in the continuum. 

By means of a power search across radial velocity data. 
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Figure 6. V/R ratio versus HJD for H/3 and Hell A4686 A. We find QPOs during all three nights. 



paper I determined that the characteristic velocities of the 
QPO signal were intermediate between the orbital regime 
and the spin regime. This is consistent with a signal ori- 
gin in the inner disc or the threading region between disc 
and magnetic curtain. However the search had no means 
to discriminate between blue- and red-shifted material and 
therefore could say little more about the QPO mechanism. 
In this paper we conduct a similar power search but in the 
line fiuxes across discrete velocity bins of the emission pro- 
files. The resulting power maps for H/3 and He ll A4686 A are 
plotted in velocity-frequency space and provided in Fig. 10, 
and Table 7 lists the QPO period on differentq nights sam- 
pled at —400 km s^^ and — 600kms~^ in the line profile. 

We find that power associated with the QPO is not 
symmetrically distributed about the rest wavelength, but 
bieised towards the blue wing of each line, as we have previ- 
ously noted in Sec. 3.3. QPO power extends from —500 to 
— 1 000 kms~^ depending on the line and night but the QPO 
does not appear to be a strong function of velocity consis- 
tent with our results from paper I. We also find, the QPO 
tending to longer frequencies with time, as we have already 
determined from paper I and the V/R ratio analysis in the 



current paper. We discuss the significance of this result in 
Sec. 4. 

We also find power on the white dwarf spin frequency 
extending to ~ 1000 km s~^, although weaker than we have 
found in the radial velocity analysis of paper I. 



4 DISCUSSION: QPO MECHANISMS 

In paper I we found that the optical QPOs observed during 
the 1996 outburst have a velocity structure that is consis- 
tent with a mechanism where dense blobs of gas orbit in the 
inner disc at a radius determined by the impact of an over- 
flowing gas stream (HL). We determined that the QPO is an 
approximately constant function of velocity ruling out mech- 
anisms involving radial or vertical oscillations in the inner 
disc flow (Carroll et al. 1985), and that the QPO tends to 
higher frequencies with time. In this paper we have shown 
that the QPO is biased towards blue-shifted material and 
we discuss this result in terms of the disc-overflow accretion 
model and the alternative beat model proposed by WKO. 
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Figure 7. Power spectra obtained from the V/R ratios of H/3 and HellA4686A on 1996 Feb 26 to 28 after the QPO has been filtered 
out. 



Table 6. The QPO period determined from the integrated flux 
of H/3 and He II A4686 A and for the continuum the three nights 
of observations. The error given is the minimum error due to the 
frequency sampling. The third column gives an estimate of the 
confidence of the identification. This estimate has been obtained 
by using a randomisation Monte-Carlo Technique over 1000 per- 
mutations. 



1996 Feb 



■PqPO (s) 



Percentage 
of Confidence 



H/3 26th 


5932 ± 2 




89 




H/3 27th 


4221 ± 1 & 


3246 ± 1 


95 & 


c 95 


H/3 28th 


3596 ± 1 




93 




Hen 26th 


6311 ± 2 




90 




Hen 27th 


4435 ± 1 & 


3228 ± 1 


95 & 


I 95 


Hen 28th 


3447 ± 1 




95 




Continuum 26th 


6750 ± 26 




70 




Continuum 27th 


4477 ± 12 




95 




Continuum 28th 


3692 ± 9 




94 





4.1 The Disc-overflow accretion model 

This model was proposed by HL on the basis that a 5000 s 
period is consistent with the Keplerian period of blobs of 
material deposited in the inner disc by the overflowing gas 
stream (Lubow & Shu 1975; Lubow 1989; Hellier 1993; Ar- 
mitage & Livio 1996, 1998), and that the X-ray hardness 



Table 7. The QPO period determined from bins across the emis- 
sion line profiles of flux of H/3 and He II A4686 A on the three 
nights. Errors are, as in Table 6, the estimates of the minimum 
error. All periods have been found to be significant to 95 per cent. 



1996 Feb Pqpo («) 

at -400 km s' 



■PqPO (s) 

at —600 km s" 



H/3 26th 
H/3 27th 
H/3 28th 
Hen 26th 
Hen 27th 
Hen 28th 



6336 ± 2 
4765 ± 1 
3766 ± 1 
6445 ± 2 
5912 ± 1 
3725 ± 1 



4645 ± 1 
3637 ± 1 
6378 ± 2 

4038 ± 1 



ratio measured as a function of QPO phase from the data 
collected by WKO is consistent with the photo-electric ab- 
sorption of soft X-rays by cool intervening gas. The optical 
counterpart to the X-ray QPOs could either be direct repro- 
cessing off the blobs or the periodic reprocessing off material 
in the outer disc as the blobs intermittently absorb the X- 
rays from the central object. 

We have determined that, in velocity, the QPO ranges 
between the expected rotational velocity of the outer disc 
and the overflow impact site (paper l). Although this distri- 
bution of QPO power is consistent with the proposed mech- 
anism, it is more difficult to reconcile the blue-shifted bias 
of power in terms of the overflow model (Fig. 10). 
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Figure 8. The QPO-filtered V/R ratios versus spin pliase obtained 
over a second spin cycle. 

It is unlikely that the optical QPOs can be the result of 
reprocessing off the disc unless there is an emission mecha- 
nism within the disc which is extremely anisotropic. Simi- 
larly direct emission from the blobs must also be anisotropic. 
In this case the cooling of shocked- or viscous-heating gas 
within the blobs could cause the anisotropy provided the 
blobs are orbiting faster than the surrounding disc material. 
It is not clear why this should be the case. 

4.2 A disc-curtain beat mechanism 

WKO proposed a QPO mechanism where the observed pe- 
riods are the beats between the spin frequency of the white 
dwarf and the Keplerian frequency of dense blobs of gas or- 
biting at the inner rim of the disc (Alpar & Shaham 1985 
a, b). Each time an accretion curtain sweeps over a blob 
we observe an increase in column density across the curtain, 
providing a cool absorbing body for the X-ray emission. The 
prediction follows that the orbital period of the inner disc is 
either ~ 320 s or ~ 380 s. We investigate whether this model 
can explain the observed bias in the QPOs across the optical 



0.15 r 




-0.15 t .............. . 

0.0 0.5 1.0 1.5 

Spin Phose 



H/3 and He II A4686 A on 1996 Feb 26 to 28. The data are repeated 



emission lines. Two schematics of the accretion flow in the 
binary are depicted in Fig. 11. We take the orbital inclina- 
tion to be consistent with 46° < i < 72° (Reinsch 1994). 
We have assumed that the magnetic axis of the white dwarf 
is misaligned with the rotational axis of the system by 45°, 
although its true inclination is unknown. 

In Sec. 3.3 we determined that similar to the X-ray 
QPOs, the optical counterpart in the emission lines are the 
result of absorption. We consider two mechanisms which 
modulate the line emission by absorption over the WKO 
beat cycle. First we consider self-absorption of line flux from 
the curtain just above each threading region. Our observa- 
tion that QPO signal does not extend to as large velocities 
as the spin signal provides some justification for this as- 
sumption. We require the absorption profile to be saturated 
such that line strength is a function of both column density 
through the curtain and the velocity gradient across the flow 
(see e.g. Horne & Marsh 1986). At spin phase equal zero, 
(pspin ~ the upper accretion curtain lies behind the white 
dwarf and material fiowing along that curtain is blue-shifted. 
Conversely material in the curtain feeding the lower pole 
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Figure 9. Power spectra for the integrated fluxes of H/3, Hell A4686A and the continuum. 



is red-shifted. We should observe a blue shifted QPO sig- 
nal when the blob sweeps through the first threading region 
every few thousand seconds, increasing the column density 
along the approaching curtain, but not a red-shifted signal 
because the second threading region is obscured by the in- 
ner accretion disc. The result is a blue bias in the QPO 
signal across the emission lines at this spin phase. However 
at <j>spin ~ 0.5 the curtain geometry has rotated by 180° 
and both curtains are equally visible. But in this configu- 
ration the velocity gradient across the line forming regions 
is small compared to the (papin — case and consequently 
the amount of absorption across the line profile is smaller. 
In this way the blue bias in the signal is conserved over the 
beat cycle. 

An equally plausible alternative is that the absorption is 
of continuum light from the accretion disc behind the white 
dwarf. As before, the column density along the accretion 
curtains is modulated on the beat cycle as a blob sweeps 
through the threading regions. This provides a kilo-second 
QPO by periodic absorption which occurs when the upper 
curtain is back-illuminated by the disc at (l>spin ~ 0. Since 
there is no back-illuminating source for the lower curtain, or 



for the upper curtain when it is red-shifted, this mechanism 
provides a natural blue bias to the QPO signal. 

The beat model explains the long-timescale QPOs from 
GK Per using the pre-existing models of QPO generation. In 
these models the driving mechanism has a timescale of a few 
hundred seconds, as observed in the rest of the dwarf nova 
class of objects. The extra ingredient for GK Per is provided 
by its properties as both a dwarf nova and an intermedi- 
ate polar, where the QPO beats with the accretion curtains 
which thread the disc onto the rapidly spinning white dwarf 
to provide the observed kilo-second periods. In the above 
discussion we have considered the QPO in terms of a blob 
mechanism but the alternative models of radially-oscillating 
acoustic waves in the inner accretion disc work equally well. 
(Okuda et al. 1992; Godon 1995). Consequently we do not 
require a new physical explanation of QPOs to explain the 
phenomenon in GK Per. 

4.3 QPOs or DNOs? 

QPOs are present in dwarf novae during quiescence and out- 
burst. However the kilo-second QPOs in GK Per have to 
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Figure 10. Power maps of frequency versus line velocity obtained from H/3 and Hell A4686A emission line fluxes of GK Per. We detect 
power at the 4000-6000 s frequency predominantly across the blue wings of these lines. We detect further power at the white dwarf spin 
frequency (~ 246 cycles/day). 



date only been found when the system is in outburst (Rein- 
sch (1994) claims a tentative kilo-second detection in opti- 
cal photometry but provides no evidence) . This behaviour is 
more typical of another class of oscillations - the dwarf nova 
oscillations (DNOs) , which occur on timescales of tens rather 
than a few hundred seconds (Robinson and Nather 1979) . A 
characteristic of DNOs is that they tend to shorter periods 
as a system approaches the peak of its outburst (Patterson 
1981), similar to what we have found in the current observa- 
tions. Given that the beat model is correct, the oscillations 
found in GK Per show characteristics of both QPOs and 



DNOs. The timescales are consistent with QPOs, whereas 
their behaviour is more comparable to DNOs. 

The timescales of DNOs suggest they are driven at the 
inner edge of an accretion disc (see Warner 1995 for a re- 
view in DNOs), and since GK Per is an unusual dwarf nova 
in that the inner disc is truncated by the white dwarf field, it 
seems plausible that a DNO mechanism in this system would 
work on a larger timescale. If we are observing DNOs, the 
implication of the period increase over the three nights for 
the beat model is that the inner disc must be orbiting at 
320 s rather than 380 s. In the latter case a decrease in the 
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Figure 11. Schematics of the accretion flow in the intermediate polar GK Per. The position of the observer corresponds to orbital phase 
zero. The two diagrams are slices through a vertical plane across the binary line of centres, separated in time by half a spin cycle. 



period of the driving DNO from the disc will result in an 
increase in the observed beat period. We do not find modu- 
lations in our data with either a 320 s or 380 s period. Ten- 
tative detections of signal at 390 ± 20 s and 410 ± 13 s have 
been claimed by Mazeh et al. (1985) from optical photome- 
try during the 1983 outburst while Patterson (1981) reports 
a detection at 380 ± 20 s but never presented the result. If 
present during the quiescent state, a 380 s oscillation would 
suggest a QPOs classification. There have been no reports 
of a 320 s period in the literature. 



5 CONCLUSIONS 

We obtain new values for the systemic velocity and the ve- 
locity semi-amplitude of the secondary in agreement with 
previous authors. We conclude that there is no evidence for 
increased heating over the inner face of the donor star dur- 
ing this stage of the outburst. We find spin modulations in 
the V/R ratios of the lines but only tentatively in their inte- 
grated fluxes or the continuum. Spin power resolved across 
the line profiles extends to velocities of 1000 km s"'^, a large 
fraction of the freefall velocity of the central object. 

The detection of kilo-second QPOs across the optical 



emission line profiles of GK Per have provided an opportu- 
nity to test the mechanism behind the unique long-timescale 
QPOs in this object, which are an order of magnitude longer 
than QPOs normally observed in disc-accreting cataclysmic 
variables. We have rejected the model of HL which considers 
the direct effects of blobs orbiting at the Keplerian frequency 
of the annulus associated with a disc-overfiow impact site. 
Our favoured models consider the long QPO period to be 
the consequence of beating between more typical timescale 
QPOs or DNOs of ~ 300-400 s with the magnetic accretion 
curtain spinning with the white dwarf. Therefore we do not 
require a new model to explain these long timescales - the 
long oscillations are merely a consequence of the magnetic 
nature of the binary. 
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